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Chimera states are dynamical patterns in networks of coupled oscillators in which regions of syn¬ 
chronous and asynchronous oscillation coexist. Although these states are typically observed in large 
ensembles of oscillators and analyzed in the continuum limit, chimeras may also occur in systems 
with finite (and small) numbers of oscillators. Focusing on networks of 2N phase oscillators that are 
organized in two groups, we find that chimera states, corresponding to attracting periodic orbits, 
appear with as few as two oscillators per group and demonstrate that for N > 2 the bifurcations 
that create them are analogous to those observed in the continuum limit. These hndings suggest 
that chimeras, which bear striking similarities to dynamical patterns in nature, are observable and 
robust in small networks that are relevant to a variety of real-world systems. 

PACS numbers: 05.45.Xt, 89.75.Kd 


I. INTRODUCTION 

Synchronization is an important feature of swarms of 
fireflies |T], pedestrians on a footbridge [5], systems of 
Josephson junctions [3], the power grid [1], oscillatory 
chemical reactions [5] and cells in the heart and brain 
[6H8]. Since the pioneering work of Winfree [9] and Ku- 
ramoto |10j . mathematical models of arrays of coupled 
oscillators have been used to gain insight into the origin 
of spontaneous synchronization in a variety of different 
settings m- In addition to uniform synchronous and 
asynchronous oscillation, many networks of oscillators 
are known to exhibit a type of partial synchronization 
known as a “chimera” state [mis]. 

Chimera states are spatio-temporal patterns in which 
regions of coherence and incoherence coexist [Mj. These 
patterns have been reported in analysis and simulation of 
coupled oscillators with a variety of network topologies 
[IIIII1H12] and appear to be robust to a variety of per¬ 
turbations [33H25] . Recently, they were also observed in 
experiments with optical [2^ , chemical [30l [31] , mechan¬ 
ical [32| and electrochemical oscillators [33|. They bear 
a strong resemblance to patterns of electrical activity in 
the human brain [53H37] . 

When they coined the term “chimera” state, Abrams 
and Strogatz defined it as a state in which an ar¬ 
ray of identical oscillators splits into domains of syn- 
chonized and desynchronized oscillation. Recently Ash- 
win et al. formalised this idea into a definition of a “weak 
chimera” [38| such that one can prove the existence of, 
and investigate the stability and bifurcation of chimera¬ 
like solutions in small networks. For our purposes, we will 


use the term “chimera state” to refer to a trajectory in 
which two or more oscillators are synchronized (in phase 
and frequency) and one or more oscillators drift in phase 
and frequency with respect to the synchronized group; 
these are weak chimeras with the additional restriction 
of phase synchronization. 

In nature oscillators often experience non-local interac¬ 
tions with other oscillators that promote synchronization. 
It is natural to model the dynamics on these networks 
using finite systems of differential equations. In mod¬ 
els like the Kuramoto model, fully synchronized states 
can be understood as stationary solutions in a rotating 
frame of reference. Unfortunately, chimera states and 
other partially synchronized solutions are more difficult 
to characterize due to fluctuations in the local degree of 
synchrony, which is measured by an order parameter. As 
a result, until recently [38l|39], little progress has been 
made in analyzing chimera states for finite networks. 

Instead, theoretical investigations of chimera states of¬ 
ten replace a finite network of oscillators (discrete) with 
an infinite network (continuum) [TH [TT] [T^ HO] ■ In the 
continuum limit, the order parameter is stationary for 
a variety of synchronized and desynchronized solutions, 
including chimera states. This makes it possible to char¬ 
acterize chimera states by solving an eigenvalue problem. 
Unfortunately, rigorous analysis of the detailed dynamics 
and stability of chimera states remains difficult, even in 
the continuum limit (it is possible in certain lower dimen¬ 
sional cases SO])- So, researchers typically discretize the 
theoretical solutions from the continuum limit and assess 
the stability using numerical simulations with large en¬ 
sembles of oscillators. The implicit assumption is that 
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these two systems, one with an infinite number of oscil¬ 
lators and the other with a finite number, should behave 
in similar ways. However, at this time there is no rigor¬ 
ous justification for why this would need to be the case. 
Although there is evidence that this holds for certain cou¬ 
pling schemes [35], there are cases where chimera states 
are stable in the continuum limit but not when the num¬ 
ber of oscillators is finite. For example, Wolfrum and 
Omel’chenko jH] study a ring of oscillators with phases 
described by 

sin(0fc-0,+a) (1) 

j=k-R 

for k = 1,2,... N where N is the number of oscillators. 
For finite N, and R proportional to N they show that 
chimera states are not attracting, but instead appear as 
chaotic transients. This means one can observe a period 
of chimera-like behavior for certain initial conditions but 
the duration of this period is finite. In large systems 
{N > 50) the observation of the collapse of a chimera 
in simulation is extremely rare, but in small networks 
{N < 25) the lifetime can be quite short. It is unknown to 
what degree these observations apply to other topologies. 

Here we explore a simple network with two groups of 
oscillators. We show that in this network, the analogues 
of chimera states need not be transients — they exist as 
stable periodic orbits with as few as two oscillators per 
group. 


II. TWO GROUPS OF N OSCILLATORS FOR 
REDUCING N 



FIG. 1. Schematic diagram showing a network of 2N oscilla¬ 
tors in two groups of N coupled with strength p within group 
and v between groups. 


We start with the system studied in [T3 |26l |42l |43] 
consisting of two groups of N phase oscillators with 
Kuramoto-Sakaguchi type coupling sin((/j -|- a) = cos(vj — 
/3) for P = 7r/2 — a. We assume an inter-group cou¬ 
pling strength oi ly = (1 — A)/2 and intra-group coupling 
^ = (1 -I- A)/2, where 0 < A < 1; figure illustrates 
such a network. Let {9i}fL^ and represent the 

phases of oscillators in groups 1 and 2 respectively. If all 
oscillators have the same natural frequency w, then their 


phases are governed by 
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dt 


dpi 

dt 
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( 2 ) 


(3) 


Note that the A and P used here are the same as those 
in Abrams et al. HD. 


A. Dynamics and bifurcations for N = oo 



FIG. 2. Bifurcations of chimeras in the continuum limit N = 
oo of (A similar figure appears in [U].) The yellow 

and green shaded regions correspond to the respective regions 
of existence for stationary and breathing chimeras. 


The dynamics of §-(§ were studied in the limit N 
oo in m- Summarizing, they use the Ott/Antonsen 
ansatz |43l [44] to derive a set of ordinary differential 
equations (ODEs) satisfied by the complex order param¬ 
eters for the two groups of oscillators (which measure the 
degree of synchrony within each group). A chimera state 
in our sense corresponds to one group being perfectly 
synchronized (so that the magnitude of its order param¬ 


eter \z\ = 


A. 

N = 


i=i' 


is equal to 1) while oscillators in 
the other group are asynchronous (so that magnitude of 
its order parameter satisfies 0 < | 2 | < 1). In this case, 
Abrams et al. derive a pair of real ODEs for the magni¬ 
tude of the order parameter of the asynchronous group 
and the difference in phases of the two order parameters. 
Analyzing these equations on varying A and P, they ob¬ 
tain the behavior shown in Fig. For small P, as A is 
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increased from zero, two chimera states (one stable and 
one a saddle) are created in a saddle-node bifurcation. 
As A increases the stable chimera undergoes a supercrit¬ 
ical Hopf bifurcation, leading to a “breathing” chimera. 
Increasing A further results in this solution colliding with 
the saddle chimera in a homoclinic bifurcation, and no 
stable chimeras remain. 


B. Dynamics for finite N 

Generating initial conditions consistent with the sta¬ 
tionary chimera in the continuum limit for varying JV 
we obtain the trajectories for order parameter x = 
shown in Fig. When N is small, it ap¬ 
pears that the dynamics oirhe order parameter in the 
desynchronized group are dominated by fluctuations, a fi¬ 
nite size effect not observed in the infinite N limit. As N 
increases, the magnitude of these fluctuations decreases 
and the dynamics approach those of the continuum limit. 
For TV > 50 the dynamics of the order parameter are vir¬ 
tually indistinguishable from the results predicted by the 
continuum theory. 


N=2 N=3 N=5 



t t t 

N=10 N=50 N=~ 



t t t 


N=2 N=3 N=5 



t t t 

N=10 N=50 N==o 



t t t 


FIG. 4. The trajectories of the order parameter in the 
asynchronous group for various values of N. Parameters; 
A = 0.4,/3 = 0.025. 


abe and Strogatz showed that arrays with JV globally- 
coupled identical phase oscillators possess TV — 3 con¬ 
stants of motion, provided that TV > 3, for almost all 
initial conditions [IS]. In other words, only three gov¬ 
erning equations are needed to describe the dynamics. 
Pikovsky and Rosenblum [46] extend this analysis to net¬ 
works with multiple groups [46]. We therefore consider 
their alternative parametrization of the system in which 
the dynamics of each group are governed by 


dt 2 ^ ^ ^ 


dt 


2pj 




1 + p 


d^4 

—A =uj+ „ 
at 2pj 


lm{Zj 




(4) 

(5) 

( 6 ) 


for j = 1,2, where 


FIG. 3. The trajectories of the order parameter in the 
asynchronous group for various values of N. Parameters; 
A = 0.1,d = 0.025. 

For parameter values above the Hopf bifurcation shown 
in Fig. [^ we obtain the results shown in Fig. [^ For 
large but finite TV we see similar fluctuations superim¬ 
posed on the oscillations of the order parameter shown 
in the TV = oo case. However, such an oscillatory state 
does not appear to be stable when TV is small, and it 
is not clear why this is the case. Therefore, at least for 
small TV, the finite size effects dominate the dynamics 
and the intuition that one can gain from the TV = oo 
case is limited. 

In order to study these phenomena further it is help¬ 
ful to reduce the dimensionality of the system. Watan- 


-T(l-h - i(l - 

^1=-^- (7) 

^ -T(l-hA)e*^P2e*‘^''72-*(l-^)e*'^Pie*‘^i7i 

Z 2 —-—-——- 2 ----— 

where 

(i) 

and the '0^ , with It = 1,.. .TV, are the constants asso¬ 
ciated with the Watanabe and Strogatz transformation 
for population j. In this system, pj measures the degree 
of synchrony in group j (but it is not equivalent to the 
order parameter in Ref. HZ]). and 'i/j are related to 
the mean phase and spread of the phases of oscillators, 
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respectively, within group j. Assume we are in a chimera 
state where pi = 1. Then 71 = 1 and 

-i(l + A)e*^e*‘^i-i(l- 

^- (9) 

Z _ -»(! + - i(l - 


1 ^ zvt 

TH if 1 P2e 

Te^ :=P2l2 = ^Z^-— 
A f _ 


pse*'"" + 


( 2 ) 


fc=i + p2e 


■, ( 2 ) 


Defining A = $1 — $2 we have 


dp2 

dt 

dA 

dt 

d'i>2 

dt 


1 - pj 
4 

1 +A 


[(1 + A)r sin (^ + /3) + (1 - A) sin (A + /?)] 


— cos /3 + r 


1 + P 2 
2 p 2 


cos + /3) 


1 - A 


-T cos (^ - A + /3) + 


1 + pj 

2 p 2 


cos (A + /3) 


1 - pj 
4p2 


[(1 + A)r cos (C + /?) + (1 - A) cos (A + ,5)] 


( 11 ) 

( 12 ) 

(13) 


where T and ^ are functions of 4^2 and p 2 - Now for a uni- 
form distribution of the tpl , i.e. ipl. ' = 2Trk/N (which 
is the appropriate choice for comparison with dynam¬ 
ics on the Ott/Antonsen manifold in the A —>■ 00 case) 
Pikovsky and Rosenblum [45] showed that 


(l-pi)(-p2e--‘--)^ 

l-(-p 2 e-*^=)^ 


As A —>• 00 , we have 72 —>■ 1, ^ 0, P —>• p 2 and (11 1-(|12|) 

decouple from (13), resulting in equations (12) from |T7| 
(after a redefinition of parameters). For finite A the 
equations are coupled. Simulating (11)-(13) with appro¬ 
priate initial conditions for varying A we obtain similar 
results to those in Figs. and (results not shown) al¬ 
though as noted above, p 2 does not correspond exactly to 
the magnitude of the order parameter plotted there. In 
all of these chimera states observed, '£'2 decreases mono- 
tonically. 


C. Dynamics and bifurcations for two groups of 
A = 4 phase oscillators 


We now consider the case A = 4 in more detail. To un¬ 
derstand the dynamics we place a Poincare section, D, in 
the flow at 4^ mod 2tt = tt and record the values of p and 
A as 4^ decreases through S. We drop subscripts for nota- 
tional convenience and seek chimera states corresponding 
to periodic solutions of (11)-(13|. These correspond to 


fixed points of the first return map for E. Following this 
fixed point as A is varied and /3 = 0.1 we use numerical 
continuation to obtain the results shown in Fig. 

As A is increased from zero a stable and unstable 
chimera are created in a saddle-node bifurcation. In¬ 
creasing A further results in the stable chimera becoming 
unstable via a Hopf bifurcation. Integrating ([^-([^ for 


values of A above the Hopf bifurcation we obtain the re¬ 
sults shown in Fig. This strongly suggests that the 
Hopf bifurcation is supercritical, and that the stable so¬ 
lution created at this bifurcation is destroyed in a global 
bifurcation as A is increased even further. This sequence 
of bifurcations is remarkably similar to the one observed 



FIG. 5. Top: the value of p (positive) and A (negative) on 
the Poincare section 'i> — n for eqns. (lll-(13l. Solid: stable, 
dashed: unstable. Bottom: period of periodic orbits. (Fur¬ 
ther continuation of the higher-period branch was not possible 
for numerical reasons.) Parameters: /3 = 0.1, A = 4. 
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p 


FIG. 6. The values of p and A on the Poincare section = tt 
for A — 0.28 to A = 0.35 (inner to outer) in steps of 0.01. 
Parameters: /3 = 0.1, N — 4. 



FIG. 7. Saddle-node and Hopf bifurcation curves for N = 4. 
(See Fig.|^) 


in the infinite N limit, and following the saddle-node and 
Hopf bifurcations shown in Fig.j^we obtain Fig.[^ which 
should be compared with Fig. 

Investigations for larger N (results not shown) suggest 
that as t 00 the curves of local bifurcations in Fig. 
smoothly deform into the corresponding curves in Fig. 
However, the stability of solutions indicated in Fig. [Sfis 
with respect to perturbations of just the variables in (11 )- 
(13), not to those in the original system (§-(§. Simula¬ 
tions of (H)-© for iV = 4 (again, not shown) suggest that 
chimera solutions shown as stable in Fig.j^do correspond 
to stable chimera states in ©-©. 


inal system ©-©. In a chimera state, one group is 
completely synchronized while the other is not, so, with¬ 
out loss of generality, we look for solutions in which the 
second group is synchronized, i.e. (pi = 4> for i = 1,... N 
and consider the phase differences tpi = 9i — </>. In this 
rotating frame of reference, the dynamics of the desyn¬ 
chronized cluster decouple from the synchronized cluster 
and the phase differences ipi satisfy 


dipi 1 -I- A 

dt 2 

1-A 

H-— 


N 


COS/3 - — ^ cos {ipi - tpj - j3) 




N 


— ^ cos (V'j +13)- cos (^/>i - /3) 


3 = t 


(14) 


for i = l,2...iV. Note that, in general, this frame of 
reference does not have a constant frequency, because (j) 
satisfies 


d(j) 1 -l- A 1 — A \ /I 

-^cos/3-^^^cos(V’j+^). (15) 

i=i 

Nonetheless, for the purpose of investigating the exis¬ 
tence of chimera states, this equation can be ignored (in 
Appendix we derive conditions under which having a 
synchronized group is stable). 
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D. Dynamics and bifurcations for two groups of 
A = 3 phase oscillators 

When A < 3, the Watanabe-Strogatz ansatz used 
above no longer applies. Instead, we consider the orig- 


FIG. 8. Top: the value of 'ipi,ip 2 on the Poincare section 
ip 3 = TT. Solid: stable, dashed: unstable. Bottom: period 
of periodic orbits. (Further continuation of the higher-period 
branch was not possible for numerical reasons.) Parameters: 
N = 3,3 = 0.1. 
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FIG. 9. Vertical axis: times between successive crossings of 
the Poincare section = n for a family of attractors corre¬ 
sponding to chimeras. Parameters: V = 3,/3 = 0.1. 


Setting iV = 3 in (141 we find stable chimera states for 
some parameter values. In these states, all ipi increase 
monotonically in time, and thus we can reduce the di¬ 
mensionality of the system by again placing a Poincare 
section in the flow at, say, '03 mod 27r = tt. The chimera 
state corresponds to a fixed point of the dynamics map¬ 


ping this section to itself under the dynamics of (14|. 


Following such a fixed point as A is varied we obtain 
Fig.i As in the TV = 4 case, a stable and unstable 
chimera are created in a sadle-node bifurcation as A is 
increased from zero, then the stable one loses stability 
through a Hopf bifurcation. 

Increasing A beyond the Hopf bifurcation and record¬ 
ing the times between successive crossings of the Poincare 
section 03 = tt we obtain Fig. which suggests that 
the Hopf bifurcation is supercritical. However, the bifur¬ 
cation associated with the appearance of the apparent 
period-5 orbit and its fate as A is increased further are 
unclear. Following the saddle-node and Hopf bifurca¬ 
tions as A and /? are varied we obtain Fig. These 
curves of local bifurcations qualitatively form the same 
arrangement as in Figs. and for = oo and = 4, 
respectively. As in the A^ = 4 case, numerical investiga¬ 
tion of the full system §-(§ for A^ = 3 suggests that the 
chimera states marked as stable in Fig. do correspond 
to stable chimeras in <§-§)■ 


III. DYNAMICS FOR TWO GROUPS OF A = 2 
PHASE OSCILLATORS 


FIG. 10. Saddle-node and Hopf bifurcation curves for A = 3 
(see Fig. [^. 


a factor of (1 -I- A)/2 and consider the specihc case of 
Eq. (14) with A = 2: 


dijj 


df ^ \ “ cos( 0 i - 02 - /3)] 


77 (A) 


-{cos(0i + /3) + cos (02 + /3)} - cos(0i - 0) 


^ = \ [cOs(/3) - COS(02 - 01 - 0)] 


77 (A) 


i{cos (02 + 0) + cos(0i -f /3)} - cos (02 - 0) 


(16) 


where 77 (A) = (1 — A)/(l -|- A). 

Note that Eq. (16) is invariant under symmetries 0i 


01 -b 2TTn and 02 —>■ 02 + 27rn where n is an integer. It is 
also invariant under the symmetry ( 0 i, 02 ) —>■ ( 02 , 0 i). 
Any solution will either be preserved by these symme¬ 
tries, or mapped to a new solution by the symmetry. 
Hence, any solution above the line 0i = 02 will have a 
symmetric solution below the line as well as translates 
shifted by 2Tm in either direction. 

We use sum/difference coordinates 


01 +02 


6 = 


01 - 02 


to write (161 in the more compact form: 


^ = sin^ 5 cos 13 — 277 (A) sin a cos <5 sin (3 
dS 

— = — sin i5[sin (3 cos 6 + ri{A) sin (/3 — a)] 


(17) 


(18) 


For general identical phase oscillator networks with 
fewer than four identical phase oscillators, chimera states 
in our sense cannot occur; see for example [38]. The sim¬ 
plest case of our system (two groups of A = 2 clusters) 
in which a chimera state can occur is therefore a network 
with 4 oscillators. For convenience, we rescale time by 


A. Invariant structures in phase space for A = 2 

For A = 2 the system has an integral of the motion 
when 0 = 0 , 

L(cr, (5) := cos i5 — 77 (A) cos cr. (19) 
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This can be seen by computing 

^ = sin /3 [sin^ <5 {cos 6 + rjiA) cos a} 

— 277 (A) ^ cos 5 sin^ cr] (20) 

and observing that for /? = 0 we have ^ = 0 (in the case 
A = 0 this reduces to the Strogatz-Watanabe constant of 
the motion [iS] L{ipi,tl; 2 ) = 2sin ^ sin ^). 

In addition, the symmetry of interchanging 
means that the line (5 = 0 (which corresponds to ipi = 
■02 = (t) is invariant for all and A, /3. Along this manifold, 
the flow is governed by 


da „ , ■ 

— = —2r](A) sm a sm / 
dt 


( 21 ) 


Two fixed points lie on this manifold: p^, for which 
01 = 02 = 0, representing the fully synchronized so¬ 
lution, and Ptt, for which 0 i = 02 = tt, representing a 
solution in which oscillators within each group are syn¬ 
chronized but the two groups are tt out of phase (an¬ 
tiphase) with each other (see Fig. 11). For 0 < /? < 7r/2, 


the domain where chimera states are typically of interest, 
the fully synchronized solution po is stable while the out 
of phase solution, is a saddle. 


Linearizing (18) about the manifold (5 = 0 gives 


^ = - [ 77 (A) sin(/3 - a) + sin,5] 6 + O(d^), (22) 


and thus this manifold may be repelling for some a and 
attracting for others, depending on the values of A and 
p. For parameters of interest, there are also two pairs 
of fixed points with 0 i 0 = 02 , one pair, pi and p{, near 


(0,0) and another pair, p 2 and P2, near (tt, tt) (see Fig. 11 
Appendix . pi and p[ are mapped to one another by 
the operation ( 0 i, 02 ) 1 —>■ ( 02 ,0i), as are p 2 andp^. For 
concreteness, we say that pi and P 2 are above the diago¬ 
nal, while p[ and P 2 are below. All of these are unstable. 
Analysis of the Jacobian at these fixed points reveals that 
both Pi and p 2 are saddles when A is small. As A in¬ 
creases, p 2 becomes an unstable spiral node. Eventually, 
at a critical value of A (see Appendix]^ pi and p 2 collide 
in a saddle node bifurcation. 

In order to explore the solutions of Eqns. (16), we in¬ 
tegrate the equations for a variety of initial conditions 
satisfying 0 < 0i,02 < 27r; see Fig. For fixed /3 we 
observe two different types of long-term attracting be¬ 
haviors. When /5 = 0, all trajectories remain on level 
curves of L that either connect pairs of neutrally sta¬ 
ble fixed points or are infinitely long periodic orbits. For 
0 < /3 <C 1, most initial conditions evolve toward the fully 
synchronized state (po). However, there is a set of ini¬ 
tial conditions that give rise to a stable chimera solution 
where 0i and 02 increase indefinitely at an asymptoti¬ 
cally linear rate, and for which 0 i(<) yf 02(0- In order 
to understand this behavior better we now examine the 
(01,02) phase plane as A is varied. 



FIG. 11. Schematic of equilibria for Eq. (I 6 I. As A in¬ 
creases, equilibria pi, p 2 , p'l, and p '2 move along the curve 
(blue dashed) in the direction of the displayed arrows. Ulti¬ 
mately, both pairs of equilibria collide and cease to exist in a 
saddle node bifurcation. 


B. Phase plane analysis for N = 2 


Figure [l3| displays the equilibrium points and some of 
the relevant manifolds of pi and p{, and the direction 
of flow for fixed /3 = 0.15 and six different values of A. 



FIG. 12. Trajectories for Eq. (161 with A = 0.1 and varying p. 
Top: For P = 0, all trajectories remain on level curves of the 
invariant L (19l. The gray shaded regions indicate neutrally 
stable chimera periodic orbits while the red line is a line of 
neutrally stable fixed points. Bottom: (b) For P = 0.15 fixed 
points are indicated by black circles. Trajectories within the 
gray band are stable chimeras. 


















































FIG. 13. Phase plots for Eq. (16l. Each panel contains a phase plot for /3 = 0.15 and for various values of ^ (a) ^ = 0.03, 
(b) A — 0.0662, (c) A = 0.1, (d) A = 0.194, (e) A = 0.25, and (f) A=0.71 Cyan arrows indicate the direction of flow in the 
plane. Fixed points are marked with black circles. Curves indicate selected stable (blue) and unstable (magenta) manifolds of 
the saddle points. 


This figure reveals a variety of different scenarios. When 
A is small, all trajectories converge to po (the fully syn¬ 
chronous state). As A increases, a global bifurcation oc¬ 
curs (see Fig. [T^b)). The nearly-vertical unstable mani¬ 
fold of Pi (magenta) merges with part of the stable man¬ 
ifold of pj (blue), where by “pj” we mean, in this case, 
p[ shifted vertically by 27r. This results in the creation 
of a periodic orbit with an infinite period. Beyond this 
bifurcation, a narrow “channel” bounded by the unstable 
manifolds of pi and Pi forms. Initial conditions within 
this channel (including those along the magenta unsta¬ 
ble manifolds) cannot approach any fixed point without 
crossing one of these stable manifolds (see Fig. [T^c)), 
therefore, the lifetimes of these trajectories are infinitely 
long. Over time, they approach a periodic orbit that rep¬ 
resents a chimera state. As A continues to increase, even¬ 
tually the nearly-vertical unstable manifold of pi merges 
with a stable manifold of the image of pi under the action 
(V'i)V' 2 ) ('tpi+27r,'ip2+27r) (see Fig.[I^d)). This causes 

a second global bifurcation in which the two manifolds 
exchange orientations (see Fig. [I^e)). Here the stable 


manifolds of pi and p[ can be traced backwards indefi¬ 
nitely. Thus there are arbitrarily long transients in the 
system (which we refer to as ‘backwards chimera states’), 
but all trajectories eventually approach po. The lifetimes 
of these transients are displayed in Fig. Eventually, 
these “channels” completely vanish when pi and P 2 dis¬ 
appear as Pi and p 2 approach each other (see Fig. [T^f)) 
and undergo a saddle node bifurcation. 


C. Bifurcations for A = 2 

The bifurcations described above appear for a range of 
/3 values. As is the case for A = 3, the parameter val¬ 
ues at which these bifurcations occur can be computed 
numerically by placing a Poincare section in the flow, 
say at ipi mod 2tt = tv. Following fixed points of the map 
from this section to itself we obtain the results in Fig. [15] 
We see that the bifurcations observed when N = 2 differ 
from those found above for A = 3,4 and in Ref. [T7] for 
A —> oo. Stable periodic orbits in the flow first appear 
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FIG. 14. Foliation of the phase plane. Left and right panels 
depict the phase plane for A = 0.25 and /3 = 0.15. The white 
region exhibits no oscillation. The colored regions contain 
chimera states with finite lifetimes. The number of cycles 
before a Hxed point is reached is indicated by the color. 


at the global bifurcation involving the unstable manifold 
of pi, as shown in Fig. [T^b). As A increases, a super¬ 
critical pitchfork bifurcation of this orbit occurs and the 
stable periodic orbit splits into a pair of stable orbits and 
an unstable orbit. Eventually, both of the stable chimera 
states disappear when the second global bifuration oc¬ 
curs, as shown in Fig.[T^d). The unstable periodic orbit 
persists for larger A. 



A 



A 


FIG. 15. Top: the value of on the Poincare section tpi = tt. 
Solid: stable, dashed: unstable. Bottom: period of periodic 
orbits. (Continuation to higher periods was not possible for 
numerical reasons.) Parameters: /3 = 0.1. 


Following the bifurcations shown in Fig. 15 we obtain 
Fig. Although the bifurcations involved in creating 
and destroying chimera states are different than for the 
N > 2 cases, surprisingly, two types of stable chimeras 
still exist in a “wedge” in parameter space with one cor¬ 
ner at the origin in {A, /3) space. Although only the re¬ 
sults of numerical calculations have been shown so far, 
in Appendix [P] we present analytical calculations of the 
left boundaries of the blue and green curves in Fig. |16[ 
(Ao,0) and (Ai,0), respectively. 



FIG. 16. Bifurcations for the A = 2 case. Black: saddle-node 
of hxed points; blue: pitchfork of chimeras; red: heteroclinic 
connection between pi and pi (with one variable increased by 
27r); green: homoclinic connection between pi and pi (with 
both variables increased by 27r). Stable chimeras exist in the 
shaded regions bounded by the green and red curves, and 
the line d = 0. The blue circle (Aq ~ 0.145898) and blue 
square (Ai « 0.285714) indicate analytical predictions for 
the bifurcation points. 


IV. DISCUSSION/CONCLUSION 


Chimera states have been observed in a variety of net¬ 
works [El [IlHIIl lISHZa Ea EH EHHS] but perhaps the 
simplest network in which they are seen is one consist¬ 
ing of two groups of N oscillators each, with all-to-all 
coupling within and between groups Ha El im isms]. 
Most analysis of chimera states has been in the contin¬ 
uum limit — where N is taken to be infinite — but here 
we have investigated the opposite limit, seeing how few 
oscillators were needed in a network in order to observe 
what could be referred to as a chimera. Surprisingly, we 
found that chimera states exist and are stable even for 
N = 2, and that the bifurcation scenarios for A > 2 are 
qualitatively the same as for N = oo. 

For A > 3 we used the Watanabe/Strogatz transfor¬ 
mation and the prior calculations of Pikovsky and Rosen- 
blum [15] to derive the three ODEs (|TT|)-(13) which de¬ 


scribe the asynchronous group of oscillators. For A = 4 
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we found a bifurcation scenario which was qualitatively 
the same as that for the N = oo case. For N = 3 we 
studied the three equations (141 for the phases differences 
between the two groups and also found a bifurcation sce¬ 
nario which was qualitatively the same as that for the 
N = oo case. The N = 2 case is described by the two 
ODEs (16) and the bifurcation scenario corresponding to 
this system is quite different from all other cases, as it 
involves global bifurcations and a pitchfork bifurcation. 


Note that the chimeras discussed here are “weak 
chimeras” in the sense of [3S]. In that paper, for the 
case N = 2 and a set of phase oscillators that reduces to 
the problem considered here for a special choice of pa¬ 
rameters, attracting weak chimeras are found even in the 
limit 1 / —>■ 0 for fixed fj,. 


For chimera states to be observable they must be ro¬ 
bust with respect to hetereogeneity. Previous investi¬ 
gations have shown that they are robust with respect to 
non-identical frequencies |241I25] (in the continuum limit) 
and removal of connections within a large but finite net¬ 
work m- The derivation of many of the equations in 
this paper has relied upon the oscillators having identi¬ 
cal intrinsic frequencies, but we have verified numerically 
that the stable chimera states found here do persist when 
the intrinsic frequencies are randomly chosen from distri¬ 
butions with sufficiently small variances, for all iV > 2 
(results not shown). 

Another question that arises is whether these results 
can be extended to networks of oscillators which are not 
described by just a phase variable. Clearly chimera states 
are possible in these systems, as some of the first obser¬ 
vations of chimera states were in networks of such oscilla¬ 
tors [HKin] (also see [TTHSP] !. However, all of these stud¬ 
ies considered large (or infinite) networks, and there is 
an open question as to whether the chimera states found 
here in small networks of phase oscillators exist in small 
networks of more general oscillators. 


Appendix A: Stability of the synchronized group 


In order to show that the manifold in which one group 
is synchronized is stable, we begin by considering consid¬ 
ering equations (§-§. Moving into a rotating frame of 
reference, we let = (j)i — 4> and = 9i — (p where 


dip 

dt 


= LU — 


1-f A 


cos/3 


1-A 


N 

'^cosiipj + (3). 


3=1 


After rescaling time by a factor of (1 -I-A)/2 and defining 
77 (A) = (1 — A)/(l -I- A), we obtain 


dipi 

dt 


d^i 

dt 


1 X ^ 

cos /3 - — ^ cos{ipi - ipj - 13) 
3=1 

N N 

cos{ipj+p) — cos('0i 

3=1 3=1 


v{A) 

N 


cos/3 


1 X ^ 

— y]] cos($, - -/3) 

3=1 


v{A) 

N 


N N 

cos^pjj + P) — cos(<i)i 
3=1 3=1 


(Al) 


Ipj - P) . 
(A2) 


Note that = 0 is an invariant manifold in equation 
(A2). We now derive conditions on ipj , A, /3 such that 
this state is stable. To achieve this, we treat ipj as exter¬ 
nal forcing functions and compute the Jacobian of equa¬ 
tion (A2). It is straightforward to show that along the 
synchronized manifold ($i = 0) the Jacobian satisfies 


Jii — 


N -I 
N 


sin(^) - sinpipj + /3) 

\ / j=i 


Aij = ^sin(^), 


For compactness, we define Z = A sin{ipj + P)- In 
matrix form, the Jacobian can be expressed as 

J = sin(/3)Ljv - ri{A)ZlN 

where /yr is the N x N identity matrix and 

fN-l -1 ... -1 \ 

-I N-1 : 


Ln = 


V -1 


-1 

-1 N-1 J 


is the Laplacian of a complete graph with N vertices. 
Making use of the fact that Lyr has a single zero eigen¬ 
value and a repeated eigenvalue N (with multiplicity 
N —1), it is straightforward to show that the eigenvalues 
of J are: 


Al = -r]{A)Z, A 2 = - sin(/3) - r]{A)Z 


The eigenvector corresponding to Ai is 1 representing 
a uniform phase shift along the synchronized manifold. 
Thus the stability transverse to the manifold is deter¬ 
mined by A 2 . Rearranging, we observe that A 2 < 0 is 
equivalent to 




sin/3 

77 (A)' 


2 


2N 
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Z depends on the phases of the oscillators in the desyn¬ 
chronized group and its value cannot be determined a 
priori. However, it is clearly bounded by —1 < Z < 1. 
So the synchronized manifold is always stable for r]{A) < 
sin/3, or 


, 1 — sm /3 

A > - 

1 -I- sin /3 

In practice, however jZj remains close to 0 in a chimera 
state, so the synchronous group is stable to a variety of 
perturbations provided is not too small. 


Appendix B: Asymptotic expressions for fixed 
points Pi and p 2 

When A ^ 1, the coordinates of fixed points pi and 
P 2 can be approximated as follows: 

Fixed point pi scales like 


equation in A and /3 corresponding to the saddle node 
curve. Unfortunately, thus far we have been unable to 
find a closed form solution for the fixed points. So, in¬ 
stead, we look for a perturbative expression when A ^ 1 
and /3 ^ 1. We assume that the solution satisfies the 
following where e <§/ 1 is a small parameter 

A= 1-e 
/3 ~ e/3i -b e ^/32 -b 
i’i ~ i’lo + 


Substituting this into the determinant of the Jacobian, 
expanding in e = (1 — A), we obtain the solution curve 


P 


125 

llfT 


47 , 
—A -b 
16 
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32 


A^. 


which agrees with the numerical results displayed in 
Fig. as A —>■ 1. 


■01 ~ sin(,5) cos(/3) 


3A -b - (3 sin^ (/3) -b l) A^ 


— - (3sin"^(/3) — 38sin^{,d) — 2l) A^ 
8 


■02 ~ 02 + sin(0) cos(/3) 


A+- (sin2(/3) + ll) A2 


(55sin‘‘(/3) - 142sin2(/3) - 33) A^ 


where 0 ° satisfies the equation 


tan(/3) 


1 — cos(02) 
2 sin(02) 


Appendix D: Bifurcations of chimera for N = 2 near 

/3 = 0 


The solutions and bifurcation curves described above 
for Eq. (16 1 can be computed analytically in the limit 
0 < /3 1. We compute using the Eq. (18) in the 

sum/difference coordinates 0 and exploit the fact that 
L defined in ( |I^ is, by ( [^ , close to an invariant for this 
limit. 

In particular, for /3 = 0 the level curve L = Lq is 
invariant and we can eliminate <5 


cos 6 = rj{A) cos ct -b To, 
sin^ J = 1 — ( 77 (A) cos cr -b To)^, 


(Dl) 


Similarly, fixed point p 2 scales like 


where 77 (A) = (1 — A)/(l-bA) meaning that 


01 ~ -TT — 2 sin(,S)A 2 — 2sin(/3) cos(/3)A 
- Q sin(0)3 - sin(/3)^ As 
02 ~ TT -b 2 sin(/3)A5 — 2 sin(/3) cos(/3) A 
+ Qsin(/3)3-sin(/3)) Ai. 

These expressions are useful for identifying initial condi¬ 
tions that lead to chimera states. 


Appendix C: Saddle node bifurcation curve for N = 2 


In order compute the saddle node bifurcation curve 
shown in Fig. |16[ we now consider the Jacobian of 
Eq. (16) near pi. At the bifurcation, the determinant 
be must equal to 0 , so in theory one can find a single 



1 - (r 7 (A) COSO- -b To)^ 


(D2) 


and 

5 = A(A, cr. To) := arccos ( 77 (A) cos a + Lq) . 

All trajectories with |To| < 1 — 77 (A) will wind around 
a but not J; they are “weak chimeras” and are periodic 
orbits with a period T{Lq) that can be found in integral 
form. 

Now consider 0 < /3 <C 1 and note that 

^=/3G(a,<5) + 0(/32) (D3) 


where 

G(cr, S) = sin^ S (cos 6 + r]{A) cos a) — 2p{A)'^ cos 6 sin^ cr. 
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Hence the change in L over one period T{Lq) on Lq is 
determined to lowest order in /3 by 

pTiLo) 

A(Lo):=/ G{a,A{A,a,Lo))dt (D4) 
Jt=o 

in the sense that L(t+T(L(t))) = L{t)+pK{L{t))+0{j3'^) 
for fixed A and small /3, for any region of L such that 
T(L) is bounded. If there is a (weak chimera) periodic 
orbit within /3 of L = Lq then A{Lq) = 0. The stability 
of a (weak chimera) periodic solution Lq is determined 
by a Floquet multiplier 1 + /3M{Lo) + 0{P^), where 

M(Lo) := ^{Lo) 

so we have (for 0 < /3 <C 1) linear stability if M{Lo) < 0 
and linear instability if M{Lq) > 0. One can show that 


A(0) = 



2r]{A){l-r]{A)^)cosa ^ 

1 - Tj{Ay cos2 cr 


(the integrand has symmetry /(ct + tt) = —1(a)) meaning 
there is a periodic solution close to the curve Lq = 0 that 
is a weak chimera for 0 < A. We compute 


M(0) = 27r 


/3VA- 1 - a\ 

[ ^ ) 


which has a unique zero for 0 < Aq at 


Ao := 


7-3^5 

2 


0.1458980337. 


(D5) 


For 0 < /3 <C 1, if 0 < A < Aq then M(0) < 0 and 
the periodic orbit is stable while for A > Aq it becomes 
unstable at a pitchfork of (weak chimera) limit cycles. 

The homoclinic connection between pi and pi (with 
both variables shifted by 2 'k ) for 0 < /3 ^ 1 can similarly 
be found by finding a value of A such that the chimera 
periodic orbit is at the level Lg = 2A/(1 + A) correspond¬ 
ing to the heteroclinic orbit of the integrable system at 
/3 = 0. On this orbit we can similarly calculate 


2A 2AV2 + 2A- 6 AVA 

TTaJ ~ VA-b A2 


and so in the limit 0 < /3 1 the bifurcation is at 

A(2A/(1 + A)) = 0, meaning that 

Ai = ^ « 0.285714. (D6) 

The heteroclinic connection between pi and p( (with one 
variable increased by 2Tr) is at A = 0 in the limit 0 < 
/3 <C 1, which is where the Lq = 0 chimera solution has 
unbounded period. The values of Aq and Ai are shown 
in Fig. 
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